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^ ' We review the spectral flow techniques for computing the index of the over- 

lap Dirac operator including results relevant for SUSY Yang-Mills theories. We 
\ describe properties of the overlap Dirac operator, and methods to implement 

it numerically. We use the results from the spectral flow to illuminate the dif- 
ficulties in numerical calculations involving domain wall and overlap fermions. 
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I. Overlap and domain wall Dirac operators 

O 

o 

O ' In these proceedings, we review some basic properties of the overlap Dirac operator 

and how its index can be computed by spectral flow techniques. One of the side results is 
that for fermions in the adjoint representation of SU{N) we find evidence for fractional 
QhI topological charge. The presentation is pedagogical with the intent of illustrating the 

^ I origin of numerical difficulties in simulating overlap and domain wall fermions. Recent 

results from our work using overlap fermions can be found in references [jl|, H, ^, ^ . 
, . The massive overlap Dirac operator derived from the overlap formalism EI] is 

>< : ^ 

^ ; DM = 2 [1 + + (1 - /^)75e(Hi(m))] (1) 

where H/^(m) is a lattice hermitian Dirac- like operator describing a single fermion species 
with a large negative mass. The mass m is a regulator parameter for the theory. In this 
work, we use the hermitian Wilson-Dirac operator H^„(m) = 75-Dwiison(— ''^), although we 
have tested other fermion actions. The mass parameter — 1 < /x < 1 is related to the 
fermion mass by 

mf = Z~'fi{l + Oia')). (2) 
The propagator for external fermions is given by 

d-^(/.) = (i-/.)-Ud-;(/.)-i1, (3) 
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i.e. it has a contact term subtracted, which makes the massless propagator chirah 
{D-i(0),75} = 0. 

A massless vector gauge theory can also be obtained from domain wall fermions [0 , 
where an extra, fifth dimension, of infinite extent is introduced. In the version of ref . , 
one can show p that the physical (light) fermions contribute logdet -Ddw to the effective 
action with the 4-d action 



-Ddw — 2 



1 +yU+ (1 -/i)75tanh ( -ylogT 



(4) 



where T is the transfer matrix in the extra dimension and Lg its size. As long as logT ^ 
we obtain in the limit as — oo 

^Dw ^ ^ [1 + /i + (1 - /i)75e(- logT)] . (5) 

This is just the massive overlap Dirac operator up to the replacement — logT. It is 

easy to see that in the limit — > 0, where is the lattice spacing in the extra dimension 
(set to 1 above), one obtains — logT = (1 + 0{as)). 

II. Some properties of the overlap Dirac operator 

In many cases it is more convenient to use the hermitian version of the overlap Dirac 
operator ([1|): 

if,(/i) = 75^ov(/i) = ^ [(1 + /i)75 + (1 - l^)e{H^)] . (6) 

The massless version satisfies, 

TO0),75} = 2i/,». (7) 

It follows that [ifQ(0),75] = 0, i.e. the eigenvectors of H'^{0) can be chosen as chiral. 
Since 

^o(/^) = (l-/i')^o(0)+/x' (8) 

this holds also for the massive case. 

The only eigenvalues of Ho{0) with chiral eigenvectors are and ±1. Each eigenvalue 
< < 1 of H^{0) is then doubly degenerate with opposite chirality eigenvectors. In 
this basis Ho{^) and Dov{fi) are block diagonal with 2x2 blocks, e.g 



, (i-;,)A2 + ^ ii-f,)\VT^\ ^_fi o\ 
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FIG. 1. Plots of the optimal rational function approximation to e(x) for various order polyno- 
mials. 



For a gauge field with topological charge Q 7^ 0, there are, in addition, \Q\ exact zero 
modes with chirality sign(Q), paired with eigenvectors of opposite chirality and eigenvalue 
1. These are also eigenvectors of Ho{fi) and Do^^ji): 

-Dov (/i) zero sector : (qi) 

(10) 

depending on the sign of Q. 

We remark that from eigenvalues/ vectors of H^{0) those of both Ho{fi) and Dov{fi) 
are easily obtained. There is no need for a non-hermitian eigenvalue/vector solver! For 



example, the Ritz algorithm [10| will do just fine 



III. Implementations of the overlap Dirac operator 

In practice, we only need the application of D{^) on a vector, D{fi)ilj, and therefore 
only the sign function applied to a vector, e(if^)'?/'. Since we need the sign function of an 
operator (a large sparse matrix) this is still a formidable task. 

Methods proposed for this computation are: 

• A Chebyshev approximation of e{x) = over some interval [6, 1] [|n|]. For small 
6 a large number of terms are needed. 

• A fractional inverse method using Gegenbauer polynomials for This has a 
poor convergence since these polynomials are not optimal in the Krylov space. 



3 



1.0 



0.5 



^ 0.0 



-0.5 



-1.0 



0.3 



0.2 



0.1 



0.0 





m 


= 0.41 




o^-^o □ m 


= 3.21 




y V + m 


= 2.23 " 




L A X m 


= 3.95 - 






g . 
• t . 


o / 







2.0 



1.5 



B 1.0 



0.5 



0.0 



2 4 6 8 

t 



FIG. 2. Results for a single 8^ instanton with radius p = 2.0 and Dirichlet boundary conditions. 

Left: spectral flow of H/,(m), center: profile of the zero modes, right: mass crossing 
value as a function of the instanton radius. 

• Use a Lanczos based method to compute based on the sequence generated for 
the computation of ^ [|13| . 



Use a rational polynomial approximation for e{x) which can then be rewritten as a 
sum over poles: 



e(x) <— X 



X 



The application of x ^ ^{Hw)'ip can be done by the simultaneous solution of the 
shifted linear systems 



(12) 



{Hi + hk)(l)k = ^, X = H^ico^ + E ^k(f)k)- 

k 



One such approximation, based on the polar decomposition W^, was introduced 



in this context by Neuberger ||16[. We use optimal rational polynomials [|T|. The 
accuracy of this approximation is shown in Fig. 

We note that in all methods listed above, one can enforce the accuracy of the 
approximation of e{x) for small x by projecting out the lowest few eigenvectors of Hy^ and 
adding their correct contribution exactly. 
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FIG. 3. On the left p{0; m) of H^(m) for quenched Wilson /? = 5.7, 5.85 and 6.0. On the right, 
the approach of p(0; m) to the continuum limit in the quenched theory at fixed masses. 



To invert D'^D for overlap fermions, we have, generically, an outer CG method (a 4-d 
Krylov space search) and an independent inner search method for e{H^)ijj - maybe CG 
again. For domain wall fermions, on the other hand, a 5-d Krylov space search method 
is used. It may pay off to try to combine inner and outer CGs for overlap fermions by 
reformulating them into a 5-d problem [^, 

IV. Index defined via the Overlap formalism 

The massless overlap Dirac operator is 

^ov = ^ [1 + 75e(H,.)] ; H,. = ^M-m) . (14) 

The index is given by Q = tr e(HL)/2. We see Q simply counts the deficit of the number 
of positive energy states of Hl- A simple method of computing Q at some fixed m is via 
the spectral flow method Consider the eigenvalue problem 

IIi(m)0fc(m) = \k{m)(j)k{m), ^ -0l("^)750fe(H ■ (15) 

An efficient way to compute Q is to compute the lowest eigenvalues of H2,(m) for m > 0. 
We can prove the number of positive and negative eigenvalues of H2,(m) for m < must 
be the same, so we slowly vary the mass m from m = while keeping track of the levels 
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FIG 4 

Zero mode size distribution in lattice units for quenched 16^ SU(2) Wilson gauge action. 

crossing zero and direction of crossings up to some m. In this way, we get the topological 
charge as a function of m. 

We note the mass m must be greater than the usual critical mass of Hiifn), otherwise 
no topology change occurs and the overlap operator does not describe a massless chiral 
fermion. This critical mass value shifts from its free field value to some positive value for 
non-zero gauge coupling. We should also choose a mass below 2 so that in the continuum 
limit there are no doubler contributions. 

In Fig. H, we show spectral flow results for a smooth background field of a single 
instanton |T^]. There is a reflection symmetry about m = 4, namely the spectrum for 



8 — m is opposite that of m. We see a mode crosses down in eigenvalue, then crosses 
up again near 2. There is a degeneracy of 3 for the modes just beyond 2. Hence, as we 
increase m we find a sequence of the index Q of 1, —3, 3, —1 and for m = 1, 3, 5 and 
8. The one-dimensional profile for the modes associated with the crossings are 

zit)=Y.(l)lin,t)Mn,t) (16) 

n 

which we compare with the continuum solution in the center panel of Fig. ^ 

^(t;c,p)= . .(L,),,,3/2i • (17) 
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FIG. 5. Spectral flow and corresponding profile for SU(2) configurations with 7^ 4If where 
la is the number of crossings in the adjoint rep. and If is the number of crossings in 
the fund. rep. There is a degeneracy of 2 in all crossings in the adjoint rep. 

Also shown is the zero crossing point m as a function of p. We see that as the size of the 
instanton decreases the crossing point (the mass) increases. 

As we turn on gauge fields, the picture of the flows complicates and we can find 
crossings throughout the mass region beyond the critical mass Since we are 

interested in the zero modes at the crossings, we compute the density of zero eigenvalues 
p(0; m) of H^(m) by fitting linearly to the integrated density 

f\{X')dX' = p{0)X + lpiX' + ... (18) 
Jo 2 

In Fig. H we show p(0; m) for quenched SU(3) lattices. We see that for m beyond the 
critical mass region, the density p(0; m) rises sharply to a peak, then drops but is never 
zero, hence the spectral gap is always closed. A similar result is also found for two flavor 
dynamical fermion backgrounds |^ (simulated with positive physical quark mass, e.g. 
not simulated in the super-critical mass region). From a size distribution, we observe 
the zero modes are on the order of one to two lattices spacings for m beyond the main 
band of crossings, namely for m in the "fiat" region of p{0;m). We find that a physical 
quantity, like the topological susceptibihty, appears constant within errors in this "fiat" 
region, indicating that these small modes make no physical contribution. 
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FIG. 6. Probability p{A) versus A for two gauge ensembles where A 



For topology to change in a gauge field evolution, we must create dislocations. 
These produce the small modes observed above which force the spectral gap of H^„(m) to 
be closed. In the right panel of Fig. |^ an empirical fit of the density to an exponential 
of the inverse lattice spacing is shown. This result implies the density is only zero in the 
continuum limit. The density of zero eigenvalues of H^{m), p(0;m), is non-zero in the 
quenched case, but decreases rapidly with decreasing coupling [Q. Very roughly, we find 
p(0; m)/(T'^/^ ~ e~^'^ . We note the gauge action can be modified to reduce p(0;m), or 
even eliminate it altogether at some fixed m 



In Fig. ^, we show the size distribution for quenched 16^ SU(2) ensembles. The 
profile from Eq. (|TB|) is used to define a size motivated by the t'Hooft zero mode in 
Eq. (0) 



Pz 



(19) 



We see there is a large number of modes about 1 to 2 lattice units in size. There is a 
corresponding secondary peak around 5 lattice units and all the distributions are bounded 
in size. This result indicates that the size distribution of zero modes does not show 
evidence for a peak at a physical scale (as suggested by some models) even after we 



remove the small modes which are most likely lattice artifacts |20]. Instead, the observed 
scaling in lattice units is suggestive of a finite volume effect. 
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FIG. 7. Spectral flow of the hermitian domain wall Dirac operator H£)\y{m) and Wilson H^{m) 
on a single instanton background. Shown are 5D extents of Ls = 8, 16, 24 and 32. For 
Lg = 32, only at a mass separation of 0.2 units from the crossing has the eigenvalue 
dropped to 10^'^. 

V. Evidence for fractional topological charge 

In a continuum background field with topological charge Q, the index of the massless 
Dirac operator in the adjoint representation is equal to 2NQ [|^, ^ . Classical instantons 
carry integer topological charge and can thus only cause condensation of an operator with 
2N Majorana fermions. Witten argued that a bilinear gluino condensate exists in SUSY 
YM. Self-dual twisted gauge field configurations, with fractional topological charge 
exist (t'Hooft), and could explain a bilinear gluino condensate. What about non-classical 
gauge field configurations? 

The adjoint representation is real =^ spectrum of is doubly degenerate: adjoint 
index can only be even. Are all even values realized, or only multiples of 2N ? To this 
end, we studied the flow in the adjoint representation on two SU(2) ensembles at the same 
fixed physical volume We do find configurations with = 41 f (number of adjoint 
and fundamental crossings), but we also find configurations with la 7^ 4Jj. An example 
is shown in Fig. |^. 

To check if this evidence for fractional topological charge is a lattice artifact, we 
plot A = J(j — 41 f in Fig. ^. Note that A takes on only even values. The probability of 
finding a certain value of A, p(A), is plotted for two ensembles in Fig. We find that 
p{A) for |A| > 2 decreases as one goes toward the continuum limit at a fixed physical 
volume. However, p(±2) does not decrease indicating that it might remain finite in the 
continuum limit. 
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VI. Main problem for Overlap and Domain Wall fermions 

The existence of small eigenvalues illustrated in Fig. |^ hampers the approximation 
accuracy and convergence properties of implementations of e(if^). Eigenvector projection 
both increases the accuracy of the approximation and decreases the condition number, 
e.g. of the inner CG. 

The existence of small eigenvalues has implications also for domain wall fermions. 
One can show that the spectrum of — log T(m) of Eq. (1) around zero is the same as 
the spectrum of Hyj{m) 0. While the small eigenvalues of — logT(m) don't appear to 
cause algorithmic problems for domain wall fermions, they can induce rather strong Lg 
dependence of physical quantities, and hence causing the need for large Lg. 

VI-1. Domain Wall and Overlap-Dirac operator spectral flow for smooth SU(2) 

As an illustration of the effects of low lying modes of if^(m), we show in Fig. |^ the 
spectral flow of the hermitian domain wall operator HDy^{m) = ^5DD^^{m) on a smooth 
4^ X 8^ single instanton background with Dirichlet boundary conditions (BC). The 
includes a parity operator. Also shown is the corresponding hermitian Wilson H^im) 
spectral flow. The "zero" DWF eigenvalue sets in slowly in Lg. We need Lg ~ 1/Amin 
where Amin is the lowest eigenvalue of Hwim) (which is similar to the lowest eigenvalue 



10 



of log(T(m))) for e(— logT(m)) ^ tanh ^ logT(m)j. In fact, the (almost) chiral zero 
mode eigenvalue Adw("^) ~ const x (1 — tanh(Amin(m)Ls/2)) for m beyond the crossing 
indicating a sensitivity to the hermitian Wilson operator eigenvalue. 

In Fig. 1^, we show a similar plot of the hermitian overlap Dirac operator Ho{m) on 
a smooth 8"^ single instanton background with Dirichlet BC. There are strict zero modes 
after the crossing, m = 0.6, 0.7, and 0.8. Also shown are the zero mode profiles for these 
masses which are quite similar and nicely follow the t'Hooft zero mode solution. 



VII. Conclusions 



The overlap Dirac operator has the same chiral symmetries as continuum fermions. 
It has exact zero modes in topologically non-trivial gauge fields. It is therefore well suited 
for a study of the interplay of topology, with its associated exact zero modes, and chiral 
symmetry breaking, determined by the density of small eigenvalues. 

The creation of dislocations necessary for change of topology causes numerical diffi- 
culties with overlap and domain wall simulations. These dislocations are purely a property 
of the gauge actions used and are not a fundamental limitation of the chiral fermion for- 
malisms. In practice, the projection technique used in the overlap simulations is vital 
to precisely control the adverse influence of the low lying zero modes of the hermitian 
Wilson operator. The same technique can be used for domain wall simulations, but it is 
more cumbersome. Further work is directed towards lowering (and possibly eliminating) 
the density of low lying zero modes. A dynamical HMC algorithm for the overlap Dirac 
operator has also been developed 0]. 
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under which SURA operates the Thomas Jefferson National Accelerator Facility. 
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